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We study QCD with 2 colour-sextet quarks as a model for walking Technicolor, using lattice 
gauge theory simulations (RHMC) at finite temperature. Our goal is to determine if the massless 
theory is QCD-like (confining, with spontaneously-broken chiral symmetry) with a slowly vary- 
ing coupling (walks) or if it is a conformal field theory. We do this by simulating the theory at 
finite temperature and observing how the coupling at the chiral-symmetry restoration temperature 
depends on the temporal extent A', of the lattice (in lattice units). If the theory is QCD-like, this 
coupling should approach zero in the large N t limit in the manner predicted by asymptotic free- 
dom. If it is conformal, this coupling should approach a finite value in this limit, i.e. the transition 
would be a bulk transition. We discuss new results at N t = 6, 8 and 12. These preliminary results 
indicate that the coupling does decrease with increasing N t , but it is unclear if this is consistent 
with asymptotic freedom. 

We also present new results at /V, = 8 for QCD with 3 colour-sextet quarks. This theory is believed 
to be conformal, and is studied for comparison with the 2-colour case. Preliminary results suggest 
that N t = 8 is still too small to access the weak coupling limit. 
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1. Introduction 



We are interested in extensions of the standard model with a composite/strongly-interacting 
Higgs sector. Technicolor theories ^ ||] - QCD-like theories with massless techni-quarks where 
the techni-pions play the role of the Higgs field, giving masses to the W and Z - show the most 
promise. It is difficult to suppress flavour-changing neutral currents in extended Technicolor while 
giving large enough masses to the fermions. Technicolor theories which are simply scaled-up 
QCD fail the precision electroweak tests m\. Walking Technicolor theories [|5| [6L [7], where the 
fermion content of the gauge theory is such that the running coupling constant evolves very slowly 
over a considerable range of length/energy-momentum scales, can potentially avoid these problems 
ig^QQl. 



QCD with 1-^ < Nf < 3jq flavours of massless colour-sextet quarks is expected to be either a 
Walking or a Conformal field theory. (First term in the /3 function is negative, second positive.) The 
Nf = 3 theory is presumably conformal. The Nf = 2 theory could be either Walking or Conformal. 
In addition, if walking, it is minimal, having just the right number of Goldstone bosons (3) to give 



masses to the W and Z. For other arguments as to why this theory might be of interest see \ \13[ |14[].) 
We simulate the Nf = 2 theory to see if it walks. We also study the Nf = 3 theory for comparison. 
We simulate these theories at finite temperature, using the evolution of the lattice bare coupling at 
the chiral transition with N t to determine if it is governed by asymptotic freedom - walking - or if it 
approaches a non-zero constant (bulk transition) - conformal. The deconfmement transition occurs 
at appreciably smaller j8 = 6/g 2 (stronger coupling), and it is unlikely to give useful information 
on QCD-like versus conformal behaviour at the N t values we use. 

We simulate the Nf = 2 theory on lattices with ./V, = 4,6,8, 12 and hope to extend this to 
larger N t . (For our earlier work on this theory see [|l5[ [l6[ [F7|] Preliminary results indicate that 
P x (N t = 12) is significantly larger than fi x (N t = 8), but by less than what the 2-loop /3-function 
would predict. We are also performing preliminary runs with /3 fixed at a value above j5 x (N t = 12) 
on lattices with fixed N s varying N t , keeping N t < N s (we thank Julius Kuti for suggesting this). 
Here we look for the transition to the chiral-symmetry restored state as N t increases. 

We simulate the Nf = 3 theory on lattices with N t = 4, 6, 8 and hope to extend this to N t = 12. 
(For earlier work, see [|T8|].) Preliminary results indicate that j5 x (N t = 8) is probably significantly 
greater than j3 x (N t = 6) which would indicate that we are not yet at weak enough coupling. 

We are using unimproved staggered quarks and a simple Wilson plaquette action for all our 
simulations. The RHMC algorithm is used to implement the required fractional powers of the 
fermion determinant. 



DeGrand et al. are studying the 2-flavour theory using improved Wilson fermions [19, 2C, 



21, 22, 23]. Most of their simulations have been at zero temperature. So far their results are 
inconclusive, although they tend to favour the conformal field theory interpretation. The Lattice 
Higgs Collaboration are also studying this theory at zero temperature using improved staggered 



quarks. 1 24, 25, 26, 27, Eq]. Their results favour a QCD-like (and hence walking) interpretation. 
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2. Simulations of lattice QCD with 2 colour-sextet quarks at finite temperature 

2.1 Simulations at N t = 6 

We have extended our simulations for quark mass m = 0.005 on a 12 3 x 6 lattice, such that the 
j3 spacing is now 0.02 in the range 6.5 < j8 < 6.7, which encompasses the chiral transition. In this 
range we have accumulated 100,000 trajectories for each j8. 

At this quark mass, it is impossible to determine the position of the chiral-symmetry restoration 
transition with any precision from the chiral condensate itself. We determine the position of this 
transition from the peak in the disconnected chiral susceptibility. This is defined by: 

Xw = f[((¥Y) 2 )-((W)) 2 } (2-1) 

where V is the spatial volume of the lattice and T is the temperature, YY is a lattice averaged 
quantity. Because we only have stochastic estimators for YY (5 per trajectory), we obtain un- 
biased estimators of (YY) 2 as the products of 2 different estimators of YY f° r tne same gauge 
configuration. 
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Figure 1: Chiral susceptibilities on a 12 3 x 6 lattice, Figure 2: Chiral susceptibilities on a 16 3 x 8 lattice, 
N f = 2. N f = 2. 

Figure [j] shows these chiral susceptibilities for our 12 3 x 6 simulations. The m = 0.005 sus- 
ceptibility shows a clear peak at j8 « 6.60. Under the assumption that the position of this peak has 
little mass dependence (which is not inconsistent with this data), we thus conclude that the position 
of the chiral transition in the chiral (m — > 0) limit is at j3 = fi x = 6.60(2) 

2.2 Simulations at N t = 8 

We have performed simulations on 16 3 x 8 lattices in the vicinity of the chiral transition. 
Quark masses m = 0.0025, m = 0.005, m = 0.01 and m = 0.02 were used, to enable extrapolation 
to the chiral limit. For 6.6 < j3 < 6.8 which spans the neighbourhood of the chiral transition we 
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performed simulations at j3 values spaced by 0.02 to enable accurate determination of the position 
of this transition. At m = 0.0025 we performed simulations of 100,000 trajectories in length for 
each (5 in this range. For each of the 3 larger masses we simulated 50,000 trajectories at each of 
the j8 values in this range. 

Again, the dependence of the chiral condensates on j8 is too smooth to accurately determine the 
position of the chiral transition, even though the decrease in the chiral condensate with decreasing 
mass does indicate that there is such a transition above which the condensate vanishes in the chiral 
limit. Thus we again turn to the chiral susceptibility, which has a peak that diverges in the chiral 
limit, to accurately determine p x . A graph of these chiral susceptibilities for each of the 4 masses 
is shown in figure ||. 

Because the j3 values used are sufficiently close together, it is possible to use Ferrenberg- 
Swendsen reweighting to determine the position of the peak for m = 0.0025 more precisely. This 
yields j3 = fS pea k = 6.690(5). Since fi peu k clearly has very little mass dependence, we deduce that 



6.69(1). Hence 



P x (N t = 8) = 6)«0.09, 



(2.2) 



consistent with the 2-loop perturbative prediction of 0.087. 

2.3 Simulations at N, = 12 

We are simulating lattice QCD with 2 flavours of light colour-sextet quarks on a 24 3 x 12 lat- 
tice, in the neighbourhood of the chiral-symmetry restoration transition. We are running at masses 
m = 0.0025, m = 0.005 and m = 0.01, to enable continuation to the chiral limit. For 6.6 < j8 < 6.9, 
we are running at a set of jSs spaced by 0.02. In this range, we have 10,000 - 25,000 trajectories at 
each j3 and m (except at j8 = 6.6, m = 0.005, where we have 62,500 trajectories), and will extend 
this to 50,000 - 100,000 trajectories. 
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Figure 3: Chiral susceptibilities on a 24 3 x 12 lat- Figure 4: Chiral susceptibilities on a 16 3 x 8 lattice, 
tice, Nf = 2. N f = 3. 
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Figure || shows the chiral susceptibilities from these runs. While the m = 0.0025 and possibly 
the m = 0.005 susceptibilities show clear peaks, the statistics are not yet adequate to make a precise 
estimate of the positions of these peaks. The m = 0.0025 graph shows a peak at j3 = fi x « 6.78. 
If this were correct, it would mean that ji x (N t = 12) — P x (N t = 8) « 0.09, compared with the 
predicted 0.12. 

2.4 Other simulations and analyses 

In addition to measuring the unsubtracted, unrenormalized chiral condensates, we have also 
used the subtraction scheme of the Lattice Higgs Collaboration, which removes much of the 
quadratic UV divergence. 

(W)«d> = (WY) ~ (mv^—iW)) (2-3) 
V dm v J mv=m 



[g5|], where my is the valence-quark mass. While this reduces the magnitude of (yy), making 
the fact that it vanishes at large j3 more obvious, it does not help significantly in the precision 
determination of f5 x at the masses we use. 

Finally, we are simulating our theory on 24 3 x N t lattices for N t < 24 at j8 = 6.9 to try to 
observe chiral-symmetry restoration as N t is increased. Here we look at both the unsubtracted and 
the subtracted chiral condensates. Our preliminary results look promising but not conclusive. 



3. Simulations of lattice QCD with 3 colour-sextet quarks at finite temperature 

We are extending our earlier simulations of QCD with 3 colour-sextet quarks to N t = 8. Our 
runs have been performed on 16 3 x 8 lattices with m = 0.005 and m = 0.01. For m = 0.005 and 
6.3 < j3 < 6.5 (the neighbourhood of the chiral transition), we currently have 20,000 to 50,000 
trajectories for each j3. Elsewhere we have 10,000 trajectories for each Q3,m). 

Figure || shows the chiral susceptibilities from these runs. While we need to complete 50,000 
trajectories at each j3 in the above range to accurately determine f5 x , it already appears that 
fhc(Nt = 8) — f}%(Nt = 6) will be appreciable. Since we believe that this theory is conformal, when 
P x (N t ) will approach a finite value as N t — > °°, this suggests that we will need a larger N t to see 
the weak-coupling limit. Even if this theory were QCD-like, the 2-loop /3 -function would predict 
Pxfflt = 8) — P%{N t = 6) « 0.0025, which would be too small for us to see, and hence inconsistent 
with our preliminary observations. 



4. Discussion and Conclusions 



We simulate the thermodynamics of QCD with 2 colour-sextet quarks on lattices with 
N t = 4,6,8,12, in the neighbourhood of the chiral transition. If chiral-symmetry restoration is 
a finite-temperature phase transition, measuring (5 X as a function of ,/V, yields the running of the 
bare lattice coupling /3 with a at a = \/{N t T x ). Asymptotic freedom would imply that fi x — > oo 
(g x — > 0) as N t — > °° and hence a — > 0. At the masses we use, the unsubtracted, unrenormalized 
chiral condensates become small at large /3s, and appear consistent with extrapolating to zero as 
m — > 0. However, they decrease too smoothly to enable an accurate determination of (3 X . Hence we 
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have extracted (5 X from the peaks in the chiral susceptibilities. We have also examined the chiral 
condensates with the subtraction scheme used by the Lattice Higgs Collaboration. While these 
subtracted condensates show more clearly that the condensates will vanish in the chiral limit at 
large enough /3, they still do not allow for an accurate determination of fi x . 

We present preliminary results indicating that f3 x increases with ./V, over the range of A^s con- 
sidered. This suggests that the theory does walk. The change in j3 x between N t = 6 and Af = 8 is 
in excellent agreement with 2-loop perturbation theory. However, the change in j5 x between N t = 8 
and N t = 12 appears to be about 25% smaller than would be predicted from the 2-loop /3-function. 
This is of concern, since this lies in the assumed weak-coupling domain (j8 > j3 x (N t = 6)). How- 
ever, we need more statistics to accurately determine p x (N t = 12), and if it remains low, we will 
need to check that this is not a finite volume effect. It is possible that the 2-loop j3 -function is inad- 
equate to describe the running of the bare coupling for unimproved staggered lattice QCD at these 
couplings. This will then require simulations at larger N t , to distinguish walking from conformal 
behaviour. 

A series of runs performed on 24 3 x N t lattices for several N t s (N t < 24) at small quark masses, 
and a fixed coupling intermediate between fi x (N t = 12) and the expected value of j5 x (N t = 24) , does 
show an increase in both the unsubtracted and subtracted chiral condensates. More work is needed 
to determine if this implies a transition to a chirally broken theory as N t increases. 

The zero temperature properties of this theory need to be studied and the results compared 
with Fodor et al and DeGrand et al. In particular we need to determine the masses of the Higgs 
and the dilaton, to see if either of these could be the Higgs-like particle observed at the LHC. For a 
recent discussion of what this entails see reference [p6[]. 

We are extending our Nf = 3 runs to N t = 8. Preliminary results indicate that there is a sub- 
stantial increase in j5 x between N t = 6 and N t = 8. This would indicate that, for this range of N t s, fi x 
does not lie completely in the weak-coupling domain, and that we will need simulations at larger 
N t . 

Other theories we plan to study include SU {2) co i our with 3 colour-adjoint (symmetric) Majo- 
rana/Weyl quarks, and SU (4) co i our with colour-antisymmetric quarks. 
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